Exclusion Statistics in a trapped two-dimensional Bose gas 
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We study the statistical mechanics of a two-dimensional gas with a repulsive delta function 
interaction, using a mean field approximation. By a direct counting of states we establish that this 
model obeys exclusion statistics and is equivalent to an ideal exclusion statistics gas. 
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The concept of fractional exclusion statistics (FES) 
proposed by Haldane in 1991 |] , has proved to be a useful 
concept to describe the statistical mechanics and thermo- 
dynamics of certain one-dimensional models. P p[ 

The basic idea of FES is that adding a number of par- 
ticles, AN, to a system, blocks Ad of the states avail- 
able for the next particle according to the linear relation 
Ad = —gAN. Intuitively this corresponds to a repulsion 
between the particles, but only very special types of in- 
teractions give rise to this type of exclusion of single par- 
ticle states. In fact all established examples of FES are 
in one dimension (or are effectively one-dimensional like 
charged particles restricted to the lowest Landau level 
by a strong magnetic field [§-£§■) Thus, the observation 
J7|,g| that, in a Thomas-Fermi approximation, a two di- 
mensional Fermi or Bose gas with short range repulsive 
interactions has the same energy and number density as 
an ideal FES gas (treated in the same approximation), 
deserves further study. It is not obvious why this kind of 
interaction should give rise to exclusion of states in the 
sense of Haldane, and it is the purpose of this note to 
provide a statistical mechanics derivation. 

We start from the two-dimensional Hamiltonian 
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that has been used as a model for atoms in Bose conden- 
sation experiments using highly assymetrical traps P,E0[ . 
The particular form of the delta function is chosen as to 
reproduce the s-wave scattering phase-shifts in three di- 
mensions, and the dimensionless coupling g is simply re- 
lated to the corresponding scattering length and the out 
of plane extension of the assymetrical trap MM ■ We shall 
assume that the temperature is sufficiently high above 
the transition temperature that the only relevant mean 
field is the density, n, and that the fluctuations are small 
enough to be ignored. We shall return to these issues at 
the end of the paper. 



Before we analyze the statistical mechanics of (Q), we 
shall give a simple thermodynamic argument for why, at 
the mean field level, we expect exclusion statistics. For 
simplicity we consider the case with a constant external 
potential V, so that the density, n, is also constant. In a 
mean field approximation, and for a fixed number of par- 
ticles, the interaction energy term in (ffl) just amounts to 
a constant shift of the energy density, 
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where kT — 1/(3 is the temperature, and At = 

\l 2-rrh 2 (3 / 'm is the thermal wavelength. We have used 
that for a free Bose gas in two dimensions the pressure 
equals the energy density and substituted the pertinent 
virial expansion as expressed in the Bernoulli numbers 
B n pT| . This expression is consistent with the system 
being an ideal FES gas in two dimensions which is known 
to have a pressure equal to the energy density, and which 
differs from a free Bose gas only by a shift \g\^ in the 
second virial coefficient fl^| . 

Let us now consider the statistical mechanics of (Q), 
and assume that the potential V is slowly varying com- 
pared with the thermal wavelength, At- We then di- 
vide the system into cells of area a 2 , where Ar < «C 
|vV/V|, and study the statistical mechanics in each cell. 
In a mean field approximation, the one-body Hamilto- 
nian in the cell £ becomes, 
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where we approximate the potential in the box with the 
constant V(fi) with ft the position of the center of the 
cell. Also n{fp) = Ng/a 2 is the mean number density 



in the cell £, with Nt 3> 1 the corresponding average 
number of particles. 

The total energy, Eg is, as usual, not simply the sum 
of the one particle energies ef of (0) , but given by 
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where the last term compensates for the double counting 
of the interaction energy. 

The number of available one particle quantum states 
in the box, dg, below some energy, e , given that there 
are already Ng particles present, follows from ph 
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with t\, in the kinetic energy. Note that g = 1 corresponds 
to free fermions, and for general g 7 this relation immedi- 
ately hints at exclusion statistics; the number of states 
in the box decreases linearly with Ng. For a Hamiltonian 
of the form (0), this is true only in two dimensions. 

We can, however, still not conclude that our system is 
identical to an ideal FES gas. Haldane's original defini- 
tion of FES was for systems with a finite dimensional sin- 
gle particle Hilbert space, but it was later generalized as 
to include ideal gases, and the corresponding distribution 
functions were derived U$M • We shall now demonstrate 
that the box Hamiltonian Hg in (m indeed describes an 
ideal FES gas. 

Because of (0) , a microstate in the box can be labelled 
by a set of integers < fci < k 2 ■ ■ ■ < k^ t , with the 
corresponding energy, 
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The exclusion properties of this system are now manifest 
since (M) implies that the pseudo energies must satisfy, 
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It should be noted that the relation (0) holds only be- 
cause the kinetic energy and the number density scale 
as the same power of the box size a. This is true for 
the present case of particles with quadratic dispersion in 
two dimensions, but also for particles in one dimension 
with linear dispersion, which is the case for anyons in the 
lowest Landau level, or equivalently, chiral particles on a 
circle with an TV 2 type interaction ||- In fact, the two 
models studied in this paper and in p] can be exactly 
mapped onto each other by identifying the ki in ([]) with 
the "pseudomomenta" introduced in J6[. 

The interaction strength, g only enters in the combi- 
nation a = gh . Recently it was shown that by taking 
the limit g —> oo, h —> with a fixed, a can be inter- 
preted as a classical exclusion statistics parameter [ [15| . 
Using the Thomas-Fermi approximation for the system 
([j]), and taking the high T limit, it is easy to show that 
the density is given by n(f) = ns{r) (l + a^nB{rj) , 
where ns(r) is the density of a non- interacting Stefan- 
Boltzmann gas in the same potential V(r). Note that all 
^-dependence is gone, and that the classical density is 
lowered because of a classical statistics effect. 

The formal proof that the exclusion property (^) corre- 
sponds to an ideal FES gas as defined in fHf] is a straight- 
forward modification of the one given in p4|| for a multi- 
species system in the fermionic representation: Going to 
a continuum description with "momenta" pi and "pseu- 
domomenta" pi defined as 
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Note that the integers ki are not proportional to the mo- 
menta; they are proportional to the kinetic energy and 
correspond to the number of states below the energy t\ . 
Of course, the real energy eigenvalues are not equally 
spaced, but we shall assume that the box is large enough 
for this effect to be negligible. 

Next introduce the quantities ki by 



h = ki + 3 E 0(k 3 - h) 



(7) 



with 9{x) being the step function, and rewrite the energy 
as a sum of "pseudo energies", q, 
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replacing the sums over k and k by integrals in the usual 
way and denoting the corresponding particle densities in 
momentum space by v(p) and p(p), respectively, one finds 
the continuum version of Eq. (M) , 



P = P + 9 / dp'p(p')d(p - p") 
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Furthermore one has to demand conservation of the num- 
ber of particles when changing variables from p to p, i.e. 
v{p)dp — p(p)dp. Combining this with Eq.(fLlJ) gives 



v{p) = 
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Inserting this into the standard expression for the bosonic 
non-equilibrium entropy, 

S=-fc(^r) 2 /"dp[i/lni/-(l + i/)]ji(l + i/)] (13) 



exactly reproduces the entropy of a ideal FES gas JU 
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from which all thermodynamics follows. 

Since we have explicitly ignored both the possibility of 
a quantum condensate and of pairing fields other than 
the density, the results of this paper (and those of Ref. 
P|) can not be used for temperatures below or in the 
vicinity of the Bose condensation transition T c . This 
is true irrespective of whether this transition is of the 
Kosterlitz-Thouless type or not pl| . Rather our results 
should be relevant in a temperature regime where the 
exclusion statistics, due to the repulsive interaction, cor- 
responds to a small correction to the ideal Bose gas. It 
is an interesting open question whether the quasi parti- 
cles above a two-dimensional bose condensate also can 
be described using exclusion statistics. To answer this 
question one would analyze the corresponding statistical 
mechanics in a more sophisticated mean field approxima- 
tion that includes effects of phase coherence and pairing 
mean fields |17| . 

We should finally point out that while our analysis was 
entirely in the context of mean field approximations, it is 
an interesting question whether or not the full quantum 
problem of a two dimensional gas with delta function pq 
also allows a description in terms of exclusion statistics in 
some range of temperatures. Since the interaction does 
not involve any dimensional parameter, this possibility is 
not excluded. 
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